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We study the spin texture of a generic helical liquid, the edge modes of a two-dimensional topo¬ 
logical insulator with broken axial spin-symmetry. By considering honeycomb and square-lattice 
realizations of topological insulators, we show that in all cases the generic behavior of a momentum- 
dependent rotation of the spin quantization axis is realized. Here we establish this mechanism also 
for disk geometries with continuous rotational symmetry. Finally, we demonstrate that the rota¬ 
tion of spin-quantization axis remains intact for arbitrary geometries, i.e ., in the absence of any 
continuous symmetry. We also calculate the dependence of this rotation on the model and mate¬ 
rial parameters. Finally we propose a spectroscopy measurement which should directly reveal the 
rotation of the spin-quantization axis of the helical edge states. 
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I. INTRODUCTION 

Over the past years, two-dimensional (2D) topolog¬ 
ical insulators (TIs) have been theoretically predicted 
and experimentally realized in various semiconductor 
heterostructuresR 2 They are two-dimensional systems 
with a band gap in the bulk material, whose electronic 
band structure can be characterized by a nontrivial topo¬ 
logical Z2 invariant.® The band gap closes at the inter¬ 
faces between 2D TIs and topologically trivial systems, 
e.g ., the vacuum, so metallic one-dimensional (ID) edge 
channels must exist at these interfaces. 

The edge states of 2D TIs have been studied in de¬ 
tail both experimentally and theoretically. Most impor¬ 
tantly, these ID liquids are helical : in a given edge chan¬ 
nel, electrons with opposite spins propagate in opposite 
directions. 3 ® This helicity, which can be observed for 
instance using the spin Hall effect,® gives rise to a num¬ 
ber of fascinating effects such as unconventional non-local 
conductances in a multi-terminal setup 7 and effective p- 
wave pairing upon inducing superconductivity.® Another 
noteworthy feature of the helical edge states is their ro¬ 
bustness to perturbations. Since electrons with opposite 
momenta are Kramers partners, weak perturbations such 
as disorder neither lead to a localization of the electron 
wave functions nor to the opening of a gap in the edge 
state spectrum, as long as the system is time-reversal 
invariant.®^ 10 This is in stark contrast to most other ID 
liquids. 

In most works, it has been implicitly assumed that 
the electron spin is a good quantum number in the un¬ 
perturbed edge channel: the very notion that opposite 
“spins” travel in opposite directions is based on this pur¬ 
ported symmetry. However, it was pointed out early on 
that the Hamiltonians governing 2D TIs do not neces¬ 
sarily possess spin symmetry. Indeed, the seminal pa¬ 
pers by Kane and Mele 3 4 already investigated the effect 
of Rashba spin-orbit coupling in graphene, which breaks 
this symmetry. In the 2D TI materials currently un¬ 
der experimental investigation, i.e., HgTe/CdTe quan¬ 


tum well^ 11 and InAs/GaSb heterostructures EH 13 bulk 
and structural inversion asymmetry both violate the ax¬ 
ial spin symmetry of the Hamiltonian. Other promising 
candidate materials such as silicene^^ and tin films— 
also exhibit intrinsic Rashba spin-orbit coupling break¬ 
ing the axial spin symmetry. Naturally, the ID edge elec¬ 
trons emerging from such 2D Hamiltonians do not have 
a well-defined spin. 

Instead, if the axial spin symmetry is broken, then edge 
eigenstates with a given momentum k are generally linear 
combinations of spin-up and spin-down eigenstates with 
respect to a fixed quantization axis. In that sense, the ef¬ 
fect of breaking the axial spin symmetry can be regarded 
as a ^-dependent rotation of the spin of the momentum 
eigenstates. For momenta close to the time-reversal in¬ 
variant Dirac point {k = 0 or k = tt), the breaking of the 
axial spin symmetry in a helical edge state can therefore 
be described using a single parameter ko, which can be 
interpreted as the momentum scale on the which spin of 
the edge electrons rotates. 

Zero-energy properties of the edge channels, e.g., the 
linear conductance at zero temperature, are not affected 
by a perturbation of the axial spin symmetry because 
they rely only on the fact that electrons with opposite 
momenta are Kramers partners. Finite-energy proper¬ 
ties, in contrast, can indeed be affected because the scat¬ 
tering amplitude between right-moving and left-moving 
states at different momenta, say if-(k ) and ^ + (fc') 
(|fc| 7^ |Jfe'|), becomes nonzero. Moreover, since the 
strength of the symmetry-breaking perturbation can vary 
in space, devices containing point contacts or tunneling 
between edge states are sensitive to a nontrivial spin-axis 
rotation. 17 

In this paper, we will therefore investigate in detail 
the spin structure of helical edge states in time-reversal 
invariant systems with broken axial spin symmetry. Us¬ 
ing analytical as well as numerical methods, we will con¬ 
sider in particular 2D TIs with translation invariance 
based on the Kane-Mele model or the Bernevig-Hughes- 
Zhang model in the presence of spin-symmetry breaking 
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terms such as Rashba spin-orbit coupling or bulk inver¬ 
sion asymmetry. Moreover, we will investigate systems 
with rotational invariance based on these models, which 
can be experimentally realized in TI disks or flakes. 

The structure of this article is as follows: In Sec. [ilj we 
introduce the notion of generic helical liquids (GHL) and 
motivate the quantities which reveal the rotation of the 
spin-quantization axis, the main effect considered in this 
paper. In Sec. m we introduce the topological insulator 
models which we are investigating throughout this work. 
In the three following sections, the spin texture of the 
helical edge states and, in particular, the rotation of the 
spin-quantization axis and its dependence on system pa¬ 
rameters, is studied for three different set-ups: (i) exact 
diagonalization for tight-binding models on nanoribbons 
(z.e., translational invariance in one direction), (ii) ana¬ 
lytical solutions for continuum models on circular disks 
(z.e., rotational invariance), and (iii) exact diagonaliza¬ 
tion for tight-binding models on di sks ( z.e., no continu¬ 
ous symmetry). Eventually, in Sec. |VII| a spectroscopic 
consideration is discussed for real space disks, before we 


summarize our findings in Sec. VII 


II. GENERIC HELICAL LIQUIDS 

Counterpropagating helical states on a given edge of a 
2D TI are related by time-reversal symmetry. For a non¬ 
interacting translation-invariant edge state, the single¬ 
particle momentum fc is a good quantum number. If we 
denote right-moving and left-moving eigenstates of the 
Hamiltonian by the operators ^(fc) and ^-(fc), these 
are linked by time-reversal, 

Tj-WT-^iP+i-k), 

T^+(fc)T- x = -V>_(-fc). (1) 

and the minus sign in the last line follows from the fact 
that the anti-unitary time-reversal operator T, when act¬ 
ing on single fermions, should satisfy T 2 = — 1. If the 
Hamiltonian is time-reversal invariant, THT~ X = H, this 
rules out many common scattering and interaction pro¬ 
cesses. For instance, single-particle potential backscat- 
tering would produce a term containing ?/’+(fc)V ; -(fc / ) + 
h.c. which is not time-reversal invariant. 

As the relations 0 are identical to those of spinful 
fermions, it is often permissible to think of the label db 
as a spin label and to replace Vv,_(^) by ?/^(fc). If 
one considers for instance elastic scattering, the observ¬ 
able properties of the system are determined by the over¬ 
lap between states with equal energy, e.g., ?/?_(—fc) and 
^+(fc). These states have zero overlap because they are 
Kramers partners. 

However, if 'ip±(k) are indeed spin eigenstates, this or¬ 
thogonality persists even for states with different mo¬ 
menta: the overlap of ^ (fc) with 'ipi(k') vanishes even for 
\k\ 7 ^ \k'\ because of the spin degree of freedom. Impor¬ 
tantly, the orthogonality of these states does not follow 


from time-reversal symmetry. There are several scenarios 
where this distinction between becomes relevant. 

For instance, electron-electron or electron-phonon in¬ 
teractions may cause inelastic scattering. In that case, 
the overlap of states ^(fc) and 'ip-(k') at different mo¬ 
menta |fc| 7 ^ \k'\ enters the physical observables. If 
these states are not spin eigenstates, their overlap can be 
nonzero even if the system is time-reversal invariant. An 
important consequence is a deviation of the conductance 
of the helical edge state from the conductance quantum 
at finite temperatures. 1 ^ Moreover, the spin structure of 
the edge state can in principle be probed by injecting 
particles with spin polarization along a given axis and 
measuring the absorption. 

Let us define a fixed but arbitrary spin quantization 
axis, and let V’j(fc) and ^j( k ) create states with momen¬ 
tum k and given spin projection along this axis. Then, 
the left- and right-moving eigenstates of the Hamiltonian 
are in general linear combinations, 


( V>tO) 
yPilk) 



( 2 ) 


where B is a momentum-dependent unitary 2 x 2 matrix. 
We can choose the spin quantization axis in such a way 
that Hq = 1 at the Dirac point k = 0. The fact that 
the pairs p0 + (fc), **/;_(—fc)] and [Vq(fc), Vh( _ k)] are both 
related by time-reversal constrains the form of B &. In 
particular, one finds that B & must be an even function of 
k: Bk = B-k- Being a unitary matrix, this means that 
up to an irrelevant phase, the low-momentum behavior 
of Bk is characterized by a single parameter fco, 


Bk 




(3) 


up to order k 2 . Physically, fco can be interpreted as 
the characteristic momentum scale for the rotation of 
the spin quantization axis. Its value directly enters the 
formula for the temperature-dependent correction to the 
otherwise quantized edge state conductance, as shown in 
Refs. 18* and Qjk Moreover, it can be measured using 
point contacts in narrow topological insulators P 7 Since 
Bk cannot be extracted directly from our numerical sim¬ 
ulations, we focus instead on fco- In order to extract 
fco, we consider the scattering between left-movers and 
right-movers. In the simplest case of a potential scat- 
terer which is short-ranged compared to the Fermi wave¬ 
length of the edge state, but long-ranged compared to 
the penetration depth of the edge state into the bulk, 
the scattering amplitude is proportional to expression 18 


B\B kl 



dx'ipl M (x,y 0 )'ip +M (x,y 0 ) (4) 


where denotes the wave function of a right 

(left) moving edge state with momentum fc, evaluted at 
position (x,y). Here, x ( y ) is the coordinate transversal 
to the edge state (along the edge state). It has further 
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been shown 18 1 that Eq. 0 behaves at long wavelengths 
as 


B\B ki 


1 - + 


k 0 2 (k\ - k\) 


(5) 


which enables us to conveniently extract fco from numeri¬ 
cal simulations. In practice, a single-particle Hamiltonian 
matrix is diagonalized for each fc individually. Within a 
numerical diagonalization of a matrix, all eigenvectors 
will be computed with an arbitrary U(l) phase attached 
to them (for a given fc, all eigenvectors have the same 
phase factor, but for different A;-values, i.e., for differ¬ 
ent diagonalizations, the phase factors will be different 
in general). In order to overcome this technical issue, 
we consider in the following the modulus of [B k B k \ f 
which is gauge-invariant and will be the central object 
considered in this paper. We define the rotation of spin- 
quantization axis (RSQA), 


JC(k u k 2 ) = 


j dxipl M (x,y 0 )'>p +M (x,yo) 
ko 2 \ki - fcl| 


(6) 


The second TI model is the Kane-Mele model® for the 
honeycomb lattice which was originally proposed to de¬ 
scribe quantum spin Hall effect in graphene. The break¬ 
ing of axial spin symmetry is realized by a Rashba spin- 
orbit term (called Rashba 1) which might be induced by 
an external electrical field. Then we also briefly discuss 
the Kane-Mele model applied to the more realistic TI- 
candidate materials silicene, germanene, and stanene - 
the silicon, germanium, and tin analogues of graphene. 
Due to the buckled lattice structure of these materials, 
an additional Rashba spin-orbit term is allowed by sym¬ 
metries (called Rashba 2 ) which also breaks the axial spin 
symmetry. 


A. Bernevig—Hughes—Zhang model 

The BHZ modeP 5 in the continuum is defined in the 
basis {|Ei, +), \Hi, +), \Ei, —), \Hi, —)} with the follow¬ 
ing matrix: 


Hbhz = fc * ( °_ k) ) (7) 


with RSQA amplitude k ^ 2 . When considering disks with 
circular shape the momentum quantum numbers k will be 
replaced by angular momentum quantum numbers j and 
an RSQA amplitude j q 2 . In the case of two-dimensional 
flakes (be., disks not possessing a rotationally invariant 
shape) we will characterize the edge states simply by the 
energies E, leading to an RSQA amplitude Sq 2 . 

The physical relevance of the RSQA © is twofold: 
first, the expression for /C(fci,fc 2 ) occurs in matrix ele¬ 
ments of local electronic Hamiltonians evaluated in the 
basis of the helical edge states. Physical quantities which 
can be related to such matrix elements, e.g., the change 
in electric or heat conductance due to impurity scatter¬ 
ing, or the local tunneling probability into the edge state, 
will depend on JC. Moreover, a numerical evaluation of 
/C allows us to determine fco, which is a quantity char¬ 
acterizing the edge state spectrum near the Dirac point. 
It can be measured using nonlocal measurements such 
as momentum-conserving tunneling between nearby edge 
states. 


III. TOPOLOGICAL INSULATOR MODELS 

In the following, we briefly introduce the topological in¬ 
sulator (TI) models which are discussed throughout this 
paper. 

The first model is the Bernevig-Hughes-Zhang (BHZ) 
modeP which was proposed to describe the TI phase of 
the HgTe/CdTe quantum wells! 11 ^ Here we consider both 
the continuum model as well as the square lattice tight- 
binding version. The breaking of the axial spin symmetry 
(i.e., the spin S z symmetry) is accomplished by adding 
a term describing the bulk inversion asymmetry (BIA). 


where + and — can be seen as effective spin up and down, 
respectively. The 2x2 sub-block associated with a single 
spin is defined as 

h{ k) = e(k)I 2 x2 + di(k)a\ (8) 


where 


e(k)=C-D(k 2 x + k 2 y ) 


yj i 


d(k)= Ak x ,-Ak y ,M - B(k 2 x + k 2 v ) 


(9a) 

(9b) 


A, B, C and D are system or material parameters, M de¬ 
scribes the band gap and cq acts on the orbital subspace 
(E,H). This model can be regularized on a two-orbital 
square lattice,^ 


e(k) = C — 2D[2 — cos (k x ) + cos (k y )\ , ( 10 ) 

d(k) = [A sin(fc cc ), —Asm(k y ), M(k)] , (11) 

M( k) = M — 2B(2 — cos(k x ) + cos(fc y )) . ( 12 ) 

Note that the lattice spacing is set to unity throughout 
the paper. The corresponding Bloch matrix Ebuz from 
Hbkz (k) can be written as^ 

Hbhz = t(A sin (k x )a x ® s z - A sin (ky)a y ( 8 ) I 2;K 2 

+M(k)oq (8) ^2x2 + £(k)l4a;4) , (13) 

where T(fc) = (e^(fe), hq(fc), e^(fe), h^(k)) T is a four- 
component spinor, s^ denotes the Pauli matrices for 
physical spin and cq for the orbital degrees of freedom. 
Already in Ref.[5]it was mentioned that a term stemming 
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from BIA will be present which gives the following con¬ 
tribution to the Bloch matrix^!, 


Hbia 


0 0 0 —A\ 

0 0 A 0 1 

0 AO 0 ’ 

V-A 00 0 / 


(14) 


Hamiltonian 

'HkM = —t Yj c \ c 0 + *^SO ^2 V ' i jc\s Z Cj 
(ij) «b» 

+ iA R ,i£c!(s x dij) z c j , (15) 

(ij) 


with a material parameter A. The BIA term is important 
here since it couples the two spin channels and leads to 
breaking of spin S z symmetry, z.e., spin is not conserved 
anymore. Another source for breaking of the axial spin 
symmetry would be Rashba spin-orbit coupling,^ which 
we do not consider here for the BHZ model. In order 
to derive the tight-binding hopping Hamiltonian in real 
space we fourier-transform (ll3|) and (14) and obtain 


ffBHZ = (c - AD) Y (eU„ + h lK) 

ia 

>) 


e t e- — hi h ■ 

1 ^za^ia L ia L ia , 


+(M-4B)^(e 

ia 

D E ( e i+xa e i<7 + e i+ya e i<7 + h \+xa h icr + h l+ya h ia) 


+ 


+B Y1 ( e \+xa e io + e l+ya e ia ~ h i+xa h ia ~ h i+ya h ia ) 
ia 

Y (~ ie I+£t h it + ie I-&t h it + ie \+xi h n ~ ie l-xl h il) 


A 

+ 2 


- h ■ — e} - h ■ 

-\-ya L ia c 'i — ya i ta 


E(4 

ia 

('■ V'.l “ '‘b.i ) + h 


) 


where q = ( 04 , 04 ) is a two-component spinor and t is 
the hopping constant, Aso the intrinsic spin-orbit cou¬ 
pling, and Ar ? i the Rashba spin-orbit coupling. More¬ 
over, h>ij = ±1 is a phase factor and the Pauli matrix 
s a denotes the physical spin, d{j is the unit vector be¬ 
tween sites i and j. We further use the following nearest- 
neighbor vectors on the honeycomb lattice, 

4l = a (~«) • 62 = 7 7) • S3 = 7 (-a) • 

In momentum space, the KM Hamiltonian reads 
Hkm = ^km 3>(fe) with the Bloch ma¬ 
trix Hkm and the four-component spinor 4>(/c) = 

(oq(fc), b^(k), Gq(fc), b±(k)) T . The Bloch matrix is given 
by 


Wkm 



7 

-g 


Xi 

\ 


-9* 

-7 

X2 





X 2 

-7 

-g 



Xt 


-g* 

7 

/ 


(16) 


Here we used the following abbreviations g, 7 , xi? and 
X 2 , which are k dependent functions: 


g = t 


^ _|_ gSSikx/2^ik y /2 g ik y 


(17) 


In order to perform the exact diagonalization on nanorib¬ 
bons, one has to partially fourier-transform the Hamil¬ 
tonian such that one ends up with HBuz(%,k y ) or 
HBUz{k x ,y)i respectively. Of course, these representa¬ 
tions can be obtained by properly combining the previous 
equations. 


7 = 2Aso 


sm 


(k y ) — 2 sin(& 7 / 2 ) cos(V3k x /2) , (18) 


Xi 



(- V3-i) + e ik y(V3-i ) 
-\-2ie^ ikx / 2 e iky / 2 


(19) 


B. Kane—Mele model 

The Kane-Mele (KM) model ^ 3 4 was originally pro¬ 
posed to describe the quantum spin Hall effect in 
graphene. Although the spin-orbit coupling in graphene 
is far too small to be observed experimentally, the KM 
model serves as toy model to study two-dimensional TIs 
on the honeycomb lattice. Moreover, the Kane-Mele pa¬ 
pers have stimulated the search for other honeycomb- 
lattice materials with possibly heavier elements in order 
to realize strong spin-orbit coupling. Examples are sil- 
icene, germanene, and stanene as discussed below. The 
tight-binding version of the KM model is governed by the 


(VS + i) + e~ iky (—VS + i) (20) 

_ ‘2j / (D — '^3ikx/‘2g—iky/‘2 

While the g term describes the semi-metallic behavior of 
graphene, the 7 term leads to the topological insulating 
phase proposed by Kane and Mele. Note that 7 was first 
advocated by Haldane 22 in order to realize quantum Hall 
effect on the honeycomb lattice. The x terms correspond 
to Rashba spin-orbit coupling and mix the spin channels 
leading to a broken axial spin-symmetry. Typically the 
Rashba term is due to an external electrical field or the 
effect of a substrate. 
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C. Silicene, Germanene, Stanene 


As mentioned above, the KM model was introduced 
to describe a TI phase in graphene, a monolayer of 
graphite. While graphene does not feature such a phase, 
climbing up in the periodic table might be a promis¬ 
ing route as spin-orbit coupling is expected to be much 
stronger. Therefore people have considered monolayers of 
silicon, germanium, and tin, called silicene , germanene , 
and stanene. In contrast to graphene, these materials 
have a buckled honeycomb structure, z.e., the two sub¬ 
lattices of the honeycomb lattice are not coplanar. This 
gives rise to an intrinsic Rashba spin-orbit interaction be¬ 
tween next-neareset neighbor sites, which is expected to 
be much stronger than the Rashba spin-orbit interaction 
coming from an external electric field or the substrate. 
The corresponding Hamiltonian is the KM Hamiltonian 
(15) with additional intrinsic Rashba interaction, 


Hr 2 = -i-A Rj2 2_j 

m) 


k'ij G (& X dij) z Cj , (21) 


where j±ij = ±1 for the A (B) sublattice and dij is the 
unit vector between the second nearest neighbor sites i 
and j. The corresponding contribution to the Bloch ma¬ 
trix is given by 


^R,2 


/ • • 

-xt ' 

V • X3 


~X3 • \ 

' X3 

•J 


( 22 ) 


with 


X3 “g A fl>2 sin (k y ) + cos(VSk x /2) sm(k y /2) 
+V3i sin(V3k x /2) cos(k y /2)) 


(23) 


Now we investigate the eigenstates which correspond to 
the helical edge states, and which are well separated from 
the bulk bands, see the green (red) dots associated with 
momentum k\ (fe) in Fig. llwhich correspond to the right 
(left) mover. Using Eq. ([6 ) we compute the RSQA which 
has the following characteristic features: (i) for k\ = k<± 
the RSQA will vanish as both eigenstates were obtained 
within the same diagonalization process and are, hence, 
by construction orthogonal; (ii) for £q = — & 2 the RSQA 
must vanish, too, in order to satisfy the Kramers the¬ 
orem; (iii) for small momenta k\ and &2 we expect the 
the RSQA to behave like /C(Aq,&2) ~ ^o 2 \^i — &||; (iv) 
if spin-mixing terms A or A R?1 / 2 , respectively, are absent 
the RSQA will identically vanish. For the BHZ model we 
show the RSQA in Fig.[3](a) and (b). It fulfills the prop¬ 
erties (i)-(iv). In addition, we show in panel (c) a fit of 
the A:i=0 and ^2 = 0 lines to the low-energy behavior 
of Eq. © which allows us to extract k Q 2 . 

The same strategy is applied for the KM model whose 
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FIG. 1. (Color online). Spectrum of the BHZ model on a 
nanoribbon with L — 24 unit cells. Parameters: A = 5, 
B = —1, M = —2, A = 0.3, C — D — 0. 


IV. ROTATION OF SPIN-QUANTIZATION 
AXIS I: TIGHT-BINDING RIBBONS 

A. Set-up 

In the following, we consider both the BHZ and the KM 
tight-binding models on a nanoribbon geometry which 
possesses open boundary conditions in one and periodic 
boundary conditions in the other direction. In order to 
achieve this one needs to partially Fourier-transform, say 
in y direction while keeping all hoppings along the x di¬ 
rection in real space. The Bloch wave functions will have 
the following form: 

V’Or, y) = k)e lky . (24) 



We obtain the projected energy spectrum in the effec¬ 
tively one-dimensional Brillouin zone as shown in Fig.jTJ 


FIG. 2. (Color online). Spectrum of the KM model on a 
nanoribbon with L — 24 unit cells. Parameters: t — 1 , Aso = 
0.2, Ar 5 i = 0.05, and Ar 5 2 — 0. 
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(d) 
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FIG. 3. (Color online) The RSQA /C(Au,fo) is shown in panels (a) and (b) for the BHZ model, in panels (d) and (e) for the 
KM model as obtained within exact diagonalization on nanoribbons. Parameters for BHZ model: A — 5, B — —1, M — —2, 
C — D — 0, L — 32 with 250 discrete k values. Parameters for KM model: t — 1, Aso = 0.2, Ar,i = 0.05, Ar j2 = 0, L = 32 
unit cells with 250 discrete k values. In panel (c) one-dimensional cuts are shown for the BHZ model corresponding to JC(k±, 0) 
(blue circles) and /C(0, fe) (green circles) as indicated by the black dashed lines in panel (b). The same is repeated in panel (f) 
for the KM model where /C(fci,7r) (blue circles) and /C(7r, fe) (green circles) are shown as indicated by the black dashed lines 
in panel (e). The red lines in panels (c) and (e) are fits to the low-energy behavior <©>• 


typical spectrum is shown in Fig. [2] where particle-hole 
symmetry is broken due to the presence of Rashba spin- 
orbit coupling. This is also reflected in a slightly asym¬ 
metric behavior of the RSQA, see Fig.[3](d) and (e). Also 
the ki = it and k^ = 7r lines shown in panel (f) violate 
the low-energy behavior property of Eq. © for larger k 
values. 


B. Results 

We extract the RSQA amplitude k q 2 by fitting 
/C(&i,&2) to the low-energy prediction (pi). Here we fol¬ 
lowed two strategies: (i) For fixed k\ (A^) we perform 
one-dimensional fits along vertical (horizontal) lines in 
/C(jfei,Jfe 2 ) (/C(fci, fc 2 )) with k q 2 as the only fitting param¬ 
eter. After doing this for various values of k\ and /c 2 we 
average over all extracted k q 2 . (ii) We perform a full 
two-dimensional fit for /C(&i, fc 2 ) and extracted k q 2 di¬ 
rectly. Both strategies agree in high accuracy. Below we 
show the behavior of k q 2 as a function of A in case of the 
BHZ model, and as a function of Aso and Ar 5 i in case 
of the KM model. It seems that the RSQA amplitude 
depends linear on |A|, 

(25) 


For the KM model the dependence seems to be linear 
in the Rashba spin-orbit coupling, but proportional to 
|^so| -1 ^ 2 f° r the intrinsic spin-orbit coupling, 

kg 2 « Ckm \ • (26) 

Vl A so| 

Cbhz and Ckm are constants. We verified these findings 
for different widths L of the nanoribbon ruling out finite 
size effects. 

We also consider the KM model with the second 



FIG. 4. RSQA k q 2 as a function of BIA parameter A in the 
BHZ model. The linear curves fit the data points very well. 
Parameters are the same as in Fig. [3] 


&o 2 ~ Cbhz |A| . 
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FIG. 5. (Color online) Parameter dependence of k q 2 in the 
KM model (a) with respect to Aso and (b) with respe ct t o 
Ar,i. The y axes have been rescaled according to Eq. (26). 
See the main text for details. Parameters used correspond to 
those in Fig.[3] 


Rashba term. That allows a d escription of the proposed 
TI phase in silicenej^^ 23 -^ germanene, and stanenef^ 
In Fig. [6] we show a representative example of RSQA in 
the presence of the second Rashba term as realized e.g ., in 
stanene (neglecting the first Rashba term). We empha¬ 
size that treating both Rashba terms hardly influences 
the results as Ar ? i is much smaller than Ar ? 2 in these 
materials? 2 ^ We further observe that the presence of Ar ? 2 
does not break particle-hole symmetry; we find, however, 
that the second Rashba term is a much weaker source 
for generating RSQA. Fig. [6] is computed for Ar ? 2 = 0.5 
which is ten times larger than the value of Ar ? i in Fig. [3] 
In silicene, germanene, and stanene, Ar ? i still is so small 
that we can safely neglect it. 

Using realistic system parameters for Silicene, Ger- 
manen^ and Stanene as obtained within ab inito calcu¬ 
lations^ we find material values for k q 2 as presented in 
Table [I] We also computed k 0 2 based on ab inito calcula¬ 
tions for a closely related band structure which was pro- 
posed to describe a possible TI phase in Na2lrOsP 6 Al¬ 
though recent experiments suggest that the true material 
is due to strong Coulomb interactions in a magnetically 
ordered Mott phased the non-interacting band struc¬ 
ture provides a bond-dependent generalization of the KM 
model. This sodium-iridate band structure crucially dif¬ 
fers, however, from the KM case as the intrinsic spin-orbit 
coupling genuinely breaks the axial spin symmetry. Con¬ 
sidering the large values we find for k q 2 makes possible 
TI phases in transition metal oxides an promising plat¬ 
form to detect large corrections to the conductivity due 
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FIG. 6. (Color online) The RSQA JC(ki, fe) for the KM model 
with second Rashba term as present in Stanene and similar 
materials. Parameters for KM model: t — 1, Aso = 0.2, 
Ar,i = 0, Ar ,2 = 0.5, L — 32 unit cells with 250 discrete k 
values. 


to rotation of the spin quantization axis. 

Here we have determined the RSQA on nanorib¬ 
bons for the non-interacting TI models. An exciting 
perspective provides the idea to study interacting TI 
modelsp^z.e., topol ogical Hubbard-like models, on cylin¬ 
der geometries? 9 ® A promising model to investigate an 
interacting generic helical liquid would be the Kane-Mele- 
Hubbard model in the presence of Rashba spin-orbit 
coupling !® 


V. ROTATION OF SPIN-QUANTIZATION AXIS 
II: CONTINUUM DISKS 

A. Set-up 

In the previous section, we worked out the effect of 
breaking the axial spin symmetry on the spin texture 
of the helical edge states. We showed that it results in a 
rotation of the spin-quantization axis as a function of mo¬ 
mentum. In most of the experimental situations relevant 



k 0 2 a 2 

kod 

Silicene 

2.4341 • 10“ 7 

2027 

Germanene 

6.8594 • 1(T 5 

121 

Stanene 

1.0282 • 10“ 4 

99 

Na2lr03 

2.6503 • 10“ 2 

6.14 


TABLE I. Values of k q 2 and ko multiplied with the a, the 
distance between unit cells, for the buckled honeycomb lattice 
materials Silicene, Germanene, and Stanene as well as for the 
transition metal oxide Na2lr03- The used system parameters 
are taken from Refs. [25] and l26l The first Rashba-term is 
absent, Ar,i = 0, and in case of Na2lr03 none of the Rashba 
terms is considered. 
































FIG. 7. Setup as explained in the main text. R denotes the 
disk radius while r the distance from the origin. The disk is 
surrounded by an infinite potential. The helical edge modes 
are located at the boundary between TI and infinite potential. 


for transport, measurement geometries like Hall bars are 
used. Here momentum is not a good quantum number 
anymore and spatial symmetries are lost. We approach 
this situation in two steps: first we try to understand 
how the breaking of axial spin symmetry is manifested 
in real space when the topological insulator sample has a 
circular symmetry, z.e., is rotationally invariant. In the 
next section, we eventually consider real space samples 
of TIs which neither possess translational nor rotational 
symmetry. 

In this section, we discuss the BHZ model in the pres¬ 
ence of BIA for a disk geometry. In order to solve this 
model analytically, we switch to the continuum descrip¬ 
tion Eq. ([7]) where the BIA term couples both spin chan¬ 
nels. We assume that the disk is surrounded by an infinite 
potential in order to guarantee that edge states cannot 
penetrate into vacuum, see Fig. [7] The edge-states are 
moving clockwise and counter-clockwise, respectively, at 
the disk edge at radius R. 

Continuum models for TIs defined on disk geome¬ 
tries were investigated previously. The BHZ model with 
a Rashba term was investigated in Ref. [31] where the 
spectrum and wave functions were computed for TI 
Aharonov-Bohm rings and disks. A TI antidot (one 
might think of a “missing disk”) described within the 
BHZ model without BIA was considered in Ref. [32j In 
Ref.[33]the KM model for silicene in a magnetic field with 
and without intrinsic spin-orbit coupling Aso was dis¬ 
cussed. Ref. [34] investigated short cylinders made from 
3D TIs. 

We start from Eq. 0 and set C = D = 0 as the 
corresponding terms only shift the total energy and do 
not influence the spin texture of the helical edge states. 
The momentum-space representation of 0 is readily 
transformed to real space by rewriting k x = —id x and 
k y = — id y , respectively. The physical spin is represented 


by the Pauli matrices s a . In real space, Eq. 0 can be 
written in a compact way as 

T~L = Hbhz + Hbia 

— 'i'A(d x o'x G) s z dydy 0 12 ^ 2 ) (27) 

+ (Af + Bd% + Bdy ) a z 0 ^-2x2 + A s y 0 a y 

In order to exploit the rotational symmetry of the system 
we introduce polar coordinates x = r cos(0), y = r sin(0). 
The constant term RLbia is not affected by the transfor¬ 
mation of coordinate system. In polar coordinates, RLbuz 
reads 


%bhz = — A[e L( t >(Jz ® Sz id r <j x 0 s z 

+ e i^,®»,k <rv0 i 2x2 ) (28) 

+ (m + Bdl -B l -1 + (J z ® I 2k2 , 


where l z = — id^ denotes the orbital angular momentum. 
We further define the total angular momentum j z as 

jz = lz - t(s z a z ). (29) 

This operator commutes with the Hamiltonian, [H,j z \ = 
0, so the eigenstates of the Hamiltonian can be chosen as 
eigenstates of j z , 


jz^j(r, <j>) = (r, (p) . (30) 

By solving this equation the angular part can be deter¬ 
mined: 


i+(r, 4 >)\ 

VjH i+(r, (f>) 

*j El -(r, 4 >) 

4>)J 


^ jEl+ {r)e^ 2 \ 


= e ij4. 


VjH 1 + (,r) e ^/ 2 
'SjE 1 -(r)e ~ l<t> / 2 




(31) 


Moreover, the periodicity of the wave function, 
\hj(r, 0) = IE/(r, 0+27r), imposes the constraint j + 1/2 E 
Z. Thus, we have completely determined the angular 
part of the wave function. 

In order to solve the radial part of the wave func¬ 
tion we apply the Hamiltonian RL to the wave function, 
(j)) = 0), where E denotes the correspond¬ 

ing energy. This yields the following set of equations: 
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(m + Bd 2 - g (j+ r y 2)2 + B^dr - E^j * jEl +(r) + Ai (~d r + J *jH 1+ (r) ~ A A> jHl _(r) = 0 (32a) 

Ai (-d r - *iE 1 + (r) + (-M - B0 2 r + B 0-W> _ B l dr _ E ^j v jHi+ (r) + A * jEl -(r) = 0 (32b) 

A V jHl + (r) +(m + Bd 2 - B ^~]/ 2 ^ + B l -d r - E^j V jEl -(r) + Ai (d r + *m-( r ) = 0 ( 32c ) 

-A + Ai (d r - ^ 3 E 1 -{r) + (-M - Bd 2 + B - B^d r - E)V jHl _(r) = 0. (32d) 


We assume that each component of the radial part of 
Tj (r, (j)) can be written as 


^jair) = VjaZji^pr) 


(33) 


where Zj is a Bessel function (either of the first kind, sec¬ 
ond kind, or a Hankel function, depending on the bound¬ 
ary conditions) and a labels the orbital and spin degrees 
of freedom, a = (E/H, ±). In the following we will often 
make use of the recursion relation for Bessel functions, 


VpZj±i(y/pr) = J -Zj(^/pr) =fd r Zj(y/pr) . (34) 

In addition we rewrite the terms proportional to B as 
follows, 

I 


di- 


3 


1 


d r = ft 


3 + 1 


d r ~ 


(35) 


Applying Eq. (35) to the Bessel function and using the 


recursion formula (|34j) yields 

z AVP r ) = -pZjiVpr) 


d 2 - E + hi, 


This relation enables us to eliminate all derivatives from 
the set of coupled equations and to simplify them. These 
equations are linear in the \Es which allows us to organize 
them in a matrix equation. Thus we can write them 
compactly as 


1 

cq 

1 

Aiy/p 

0 

-A 

) 


f'&jE 1 A 

-Ai^/p 

—M + Bp — E 

A 

0 



^ jHi + 

0 

A 

M — Bp — E 

Ai^Jp 



%Ei- 

\ -A 

0 

-Ai^p 

-M + Bp- 

e) 


UiH,-/ 


= H9 = 0. 


(36) 


To find nontrivial solutions, we require vanishing deter¬ 
minant of H, 

[{M - Bp) 2 - E 2 + A 2 + A 2 p ] 2 - 4A 2 ^ 2 p = 0 . (37) 

Solving this polynomial equation with respect to p, we 
obtain four roots p n , n = 1,2 ,3,4 (the roots become 
doubly degenerate if A = 0 since Ti is block-diagonal in 
this case). It is possible to find closed expressions for 
the p n for arbitrary parameters A , £?, M, A, and the 
energy E. Since these expressions are very long, we omit 
them here. Nonetheless we always work with the exact 
expressions for p n . 

The choice of Bessel function is fixed by the require¬ 
ments not to diverge at the origin r = 0 , thus we 
work with Bessel functions of the first kind, denoted as 
Jj(\Jp r )• Now we can write the eigenstates in the pres¬ 
ence of finite A as 


^ j (a = 


(38) 


n= 1 


( ■ (M-Bp n ) 2 -E 2 +A 2 -A 2 p n j / r— \ i(j)/2 \ 
1 2A y /p^(M-Bp n -E) J j+\WEnijC 

Jj-\ ( v ^ r ) e_ * 0/2 

(M —.Bpn) 2 —g 2 —A 2 +/l 2 p n t / rr— \ 2 

2A (M-Bp n -E) J j-±WP n ' 

( M ~ B P") 2 -- E2 + A2 + A2 Pn J, + (^>) e ^/2 


2AAy/p^ 


The paramters A, B , M, and A are chosen such that we 
stay in the topological phase. Now we force the energy 
E to be smaller than the band gap which will allow us 
to compute the helical edge states. In this parameter 
regime we find the four roots p n to be complex. We 
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impose boundary conditions such that the wave function 0. Inserting this condition for the wave functions into 
vanishes at the boundary of the disk, Tj (r = R, </>) = the previous equations leads to the requirement that the 

coefficient matrix has a vanishing determinant: 


det 


Jj-l/ 2 (\/Pl R ) 
ViJj-i/2{VPi R ) 
V £l Jj +1 /2 ( \/Pl R ) 


(2Jj+l/ 2 {^/P 2 R) 

Jj-i/ 2 (VP2 R ) 

V2Jj^l/2(y/P2 R ) 

&Jj+l/2(\/P2 R ) 


where we introduced, 

(M - Bp n ) 2 - E 2 + A 2 - A 2 p n 


Cn = 


y/Pn{M - Bp n - E ) 


(40) 


, _ (M - B Pn ) 2 - E 2 + A 2 + A 2 Pn 

VpJ 3 ( j 


T]n 


(M - Bp n ) 2 - E 2 - A 2 + A 2 p n 
(M - Bp n - E) 


(42) 


By inserting the parameters A , L>, M, A, j, and R into 
the determinant, Eq. (39) only has a solution for two 
different values of the energy, -\-E and —E, to be de¬ 
termined. Once we have found the energy of the edge 
state, inserting all parameters and the energy E into the 
wavefunction renders the corresponding eigenstates. The 
larger the radius R is the more energy levels with higher 
values of total angular momentum j fit into the band 
gap. 

We compute energies and eigenstates for various disk 
sizes R leading to the spectral plots Fig. [8](a). Note that 
the energy spectrum is doubly degenerate as expected 
for a helical liquid (Kramers theorem). In Fig. [8](a), en¬ 
ergy levels for the clockwise moving edge state are labeled 
from j = 13/2 (level with lowest energy) to j = —13/2 
(level with highest energy). For the degenerate counter¬ 
clockwise moving edge state the j labels have opposite 
sign (note that for the sake of clarity, the legend has 
been omitted for these levels). All energy levels are fully 
classified by energy and total angular momentum. At 
sufficiently large R the energy levels tend to be equally 
spaced. When R 0, there are no energy levels which 
remain in the gap. We need a minimal disk radius to ob¬ 
tain edge states. This behavior seems intuitively correct 
as helical edge states possess a finite penetration length; 
very short radii allow the clockwise and counter-clockwise 
moving modes at “opposite edges” of the disk to overlap 
and gap out. 

We computed the energy spectrum for different values 
of A. Larger values of A result in an energy spectrum 
with smaller level spacing; the energy levels also appear 


CsJj+l/2(y/P3 R ) 
J j —1 /2 ( \/P3 R ) 
mJj-l/2(y/P3 R ) 
£,zJj+l/ 2 (yfPZ R ) 


CaJ j+i/ 2 (\/pi R ) ^ 
Jj — 1/ 2 (\PP4:R) 

mJj-i/ 2 (Vpi R ) 

^4 Jj 4*1 /2 ( , \/Pa R ') J 


(39) 


r 


now for smaller radii. Note that the particle hole sym¬ 
metry remains intact at all times (independent of the 
choices for R and A). As a consistency check, we recover 
for A —>> 0 the results for the disk without BIA (z.e., for 
two uncoupled Chern insulators with opposite chirality). 

Now we will discuss the properties of the eigenstates. 
The radial density of the wave function |\Jt 7 -(r, <j) = 0)| 2 
can be obtained from Eq. (39) where (j) can be set to 0 due 


to rotational symmetry. The values of A n are determined 
from the normalization condition. In addition, they need 
to fulfill that Tj(r = R, </>) = 0. The resulting total radial 
density is shown in Fig. [8jb). In panels (c) and (d) we 
have shown the spin-resolved radial densities. The shape 
of the total radial density is determined by the one of the 
majority spin. The contribution of the minority spin is 
smaller, but it penetrates more deeply in the topological 
insulator disk and the structure is very different com¬ 
pared to the majority spin. But the resulting broadening 
of the total density is very small. Note that we reproduce 
the behavior of the minority spin (Fig.[8](d)) for the tight 
binding disks discussed in the following section and it is, 
hence, not an artifact of the continuum approximation. 


B. Results 


Once we have determined the eigenstates we are able 
to compute RSQA for right and left mover with different 
angular momentum quantum number j, see Fig. [9j The 
shape of the overlap has a striking similarity to Eq. (J6| 
derived for the cylinder geometry. By analogy, we pro¬ 
pose the following expression for RSQA, 


£( 31 , 32 ) = 


j drip'll (r, (f> 0 )ip+,j 1 (r,(l)o) 

•-2 I *2 -2 1 

Jo \3i J 2 \ ? 


(43) 


where 0o) is the state at total angular momen¬ 

tum ji moving clockwise and ' 0 -,j 2 ( r 5 M i s state at 
total angular momentum moving counterclockwise. 

We obtain again the two diagonals j\ = 22 and j% = 
—j' 2 , respectively, where the RSQA vanishes. When 
ji = — j' 2 , the edge states are forming Kramers pairs, 
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FIG. 8. (Color online) (a) Energy spectrum as a function of the disk size R. (b) Total radial density |Tj(r, (j) — 0)| 2 as a 
function of the radius r for a disk of size R — 30. In (c) and (d) the radial density is again shown for the different spin channels 
separately. Parameters used: A — 5, B — —1, M — —2, A = 0.5. Note the different scale in (d) compared to (c). 


they have the same energy, but opposite total angular 
momentum. When j\ = j2, the eigenstates are orthogo¬ 
nal by construction. 

As in the previous section, we can extract the RSQA 
amplitude j q 2 for various values of A. While the curves 
jo“ 2 (A) are different for different disk sizes R , they agree 
when multiplied by R 2 (and are R independent), see 
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FIG. 9. A representative example for RSQA /C(ji,J 2 ) as a 
function of total angular momenta ji and j ‘2 of the continuum 
disk is shown. Parameters used: A = 5, B — — 1, M — —2, 
A = 0.5, and disk size R — 50. 



FIG. 10. Dependence of Jq 2 R 2 on A for three different disk 
sizes R — 30, 50, and 70. In addition, the result of Fig.[4]is 
again shown as dashed lines to emphasize the close agreement 
between the two different approaches. Parameters used: A = 
5, B = -1, and M = -2. 

Fig. [To] Note that R 2 j$ 2 has the same dimensionality as 
k q 2 . Fig. [To] suggests that the dependence is proportional 
to |A| 3 / 2 , 


jo 2 R 2 « C B hz|A| 3 / 2 (44) 

where Cbhz is a constant. The A dependence is different 
compared to the results for the tight-binding nanoribbons 
(shown as a dashed line in Fig. [To]) . The reader may no¬ 
tice that there are several drastic differences which might 
be the source for this mismatch: (i) disk vs. nanoribbon 
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FIG. 11. (a) Energy levels and (b) the corresponding density 
of states for a KM tight-binding disk. Parameters used: t — 1, 
Aso — 0.2, Ar,i — 0.1, Ar ,2 — 0 and N = 20 x 20 = 400 sites. 



FIG. 13. Example for a honeycomb lattice disk with hexag¬ 
onal shape. We define an effective radius r as indicated by 
the different colors. The set A of sites which is highlighted in 
yellow is used to compute the RSQA as discussed in the main 
text. The same sites are used for the spectroscopy discussion. 


and, more importantly, (ii) exact tight-binding model vs. 
continuum model. In fact, it is quite surprising that the 
actual values of j$ 2 R 2 are very close to the values of k^ 2 
shown in Fig.[4j 


VI. ROTATION OF SPIN-QUANTIZATION 
AXIS III: TIGHT-BINDING DISKS 

A. Set-up 

In this section we consider the most general situation 
where neither translational nor rotational symmetry is 



FIG. 12. Representative example for densities (red dots) of a 
single edge state for the tight binding “disk” . The dot radius 
is proportional to the density. The majority spin is shown 
in panel (a) and the minority spin in (b). Parameters used: 
t — 1 , Aso — 0.2, Ar,i — 0.1, and Ar ,2 — 0. Here the disks 
possess a rectangular shape (as sketched) and N — 20 x 40 = 
800 lattice sites. Note that density for the minority spin is 20 
times amplified compared to the majority spin. 


present. We consider the tight-binding version of the KM 
model in real space, z.e., we consider honeycomb lattice 
disks (or flakes). For convenience, we choose disks which 
have the shape of a hexagon. At the end of this section, 
we also consider the rectangular shaped TI samples as 
shown in Fig.[T2] in order to demonstrate that our analy¬ 
sis does not rely on the shape of the samples. Using exact 
diagonalization of the tight-binding Hamiltonian we can 
directly access energies and eigenstates. Energy spec¬ 
trum and the corresponding density of states is shown in 
Fig. 11 The density of states is particle-hole asymmetric 
due to the presence of the first Rashba term. 

Spectrum and density of states are characterized by 
supporting only double degenerate levels in the bulk gap 
which are (almost) equally spaced. We further observe 
the remainder of the van Hove singularity at E/t = ±1 
present in the pure nearest-neighbor tight-binding model 
on the honeycomb lattice. 

If we select one of the energy levels inside the bulk gap, 
we have direct access to the eigenfunction and its density 
distributions on the lattice. We are thus able to compute 
the density on the site (z, j) via 


Pa(i, j) = |Vv(m)| 2 (45) 

which is shown in Fig. [12] for a representative example. 
As expected we observe a mixing of the spin channels in 
the presence of Rashba spin-orbit coupling. Note that 
the contribution from the minority spin is much smaller 
than for the majority spin (for our choice of parameters 
by a factor 20). The same has been observed for the 
continuum disks discussed in the previous section. 

While the wave function is mainly localized at the 
outermost sites (keep in mind that Aso = 0.2 is quite 
large), we observe again that the minority spin pene¬ 
trates slightly deeper into the bulk (the same has been 
discussed for the continuum disks, see previous section). 

In the following, we consider honeycomb lattice disks 
with a hexagonal shape as shown in Fig.[l3] While rota- 
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FIG. 14. Radial densities of a helical edges states computed 
on the yellow path A of a tight-binding disk (see Fig. 13). 
(a) Majority spin and (b) minority spin (note the different 
scale for the y axis). Parameters used: t — 1 , Aso = 0.2£, 
Ar ? i = O.lt, and Ar j2 — 0. 



FIG. 15. Representative example for RSQA JC(Ei,E 2 ) as a 
function of energie levels E\ and E 2 of the tight-binding disk. 
Parameters used: t — 1 , Aso = 0.2 and Ar ; i = 0.05, Ar , 2 = 0, 
and with disk size R— 15 corresponding to 1350 lattice sites. 


compute the RSQA. In principle, we would like to apply 
formula (43) but for the tight-binding disks the angu¬ 
lar momentum quantum number is not well-defined due 
to the lack of continuous rotational symmetry. Further¬ 
more, if we considered another disk shape (e.g., a rect¬ 
angle or a Hall bar) an approximate angular momentum 
number is even less justified. Nonetheless, as shown in 
Fig.[lTJ each energy level inside the bulk gap is doubly 
degenerate, one corresponding to the right moving and 
the other to the left moving edge mode. Therefore we 
could generalize the RSQA by labeling the levels by the 
energy E for right and left moving modes. By analogy, 
we therefore define the RSQA for tight-binding disks, 


tional symmetry is lost, this geometry comes closest to 
a circular shape. We can then define concentric “rings” 
around the center which have a hexagonal shape. For 
example, in Fig. [l3]the sites with the same color coding 
are considered to be on the same ring; they have approx¬ 
imately the same distance from the center of the disk. In 
analogy to the continuum disks, we chose 0 = 0 and con¬ 
sider only those sites which are on the yellow highlighted 
path in Fig. [l3j A few examples for the radial densities 
of helical edge states are shown in Fig. [l4j 

Due to the discreteness of the lattice the wave func¬ 
tions do not appear to be smooth (e.g., as a function of 
radius), but their qualitative behavior is very similar to 
what was found for the continuum disks. The majority 
spin is strongly localized at the sample edge and its in¬ 
tensity dominates the wave function. As observed before, 
the radial density of the minority spin penetrates slightly 
deeper into the bulk. 


B. Results 

Having realized how to obtain the radial part of the 
wave function for the tight-binding disks we are able to 


1C{E u E 2 ) = 


(fij)e A 


(46) 


2 


El - E 


: 2 | 

2 5 


where A is a set of lattice sites indicated by the yel¬ 
low region shown in Fig.[l3| and i/j± J E{hj) is the wave 
function moving in the (counter-)clockwise direction at 
energy E, evaluated on the site (z, j). A representa¬ 
tive example of RSQA for tight-binding disks is shown 
in Fig. [15] Now we can again perform similar fitting 
procedures as explained in the previous sections. We 
extract the RSQA amplitude for different parame¬ 
ters and different effective radii. It turns out that the 
situation is comparable to the continuum disks and we 
should multiply the RSQA amplitude with R 2 . In addi¬ 
tion, we assume that the parameter dependence should 
be the same as in Eq. (26). Therefore we have rescaled 


the extracted data accordingly, see Fig.[16j We find very 
good agreement of the parameter dependence compared 
to the KM model on a nanoribbon. For illustration we 
have shown one of the curves from Fig. [5] as the dashed 
line in Fig.[l6] Note that we have to rescale the data 
from the tight-binding ribbons by a factor of three be¬ 
cause we assumed translation vectors a ? on the ribbon 
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VII. SPECTROSCOPY ASPECT OF TI DISKS 


In the two previous sections we found indication of how 
the RSQA can be extracted in real space. In rotation- 
ally invariant systems we find the direct analogy to the 
nanoribbons where translational symmetry is conserved. 
The consideration of the tight-binding disks revealed, 
however, that the relevant information to compute RSQA 
can be locally obtained. For the tight-binding disks we 
showed that it is sufficient to sum up the overlap contri¬ 
butions from different sites along the yellow highlighted 
path in Fig.[l3j Since we are dealing here with edge states 
naively one might think that they are perfectly localized 
at the outermost lattice site of the disk. The numer¬ 
ical results show, however, that the helical edge states 
are typically localized on a few outer sites, although the 
dominant weight is indeed located at the edge site. In 
an idealized case it would be sufficient to measure the 
wave functions for different energies on a single site to 
obtain the RSQA directly. Using spectroscopic energy- 
resolved measurements that might be even accessible in 
experiments. Therefore we investigate now the penetra¬ 
tion depth of the edge states in the context of RSQA. 
We want to understand how many sites away from the 
edge have to be measured to obtain RSQA. Practically, 
we simply modify Eq. (46) in the following way: 


K^ n \E u E 2 ) = 


j) V’+.Ei (i,j) 

(i,j)e An 


(47) 


FIG. 16. Dependence of £q 2 R 2 on (a) the intrinsic spin-orbit 
coupling Aso and (b) the Rashba spin-orbit coupling for tight- 
binding disks (t = 1) with different disk sizes R = 10, R = 
12, and R = 15 corresponding to N = 600, 864, and 1350 
lattice sites, respectively. The dashed line in panel (a) and 
(b) correspond to one of the curves shown in Fig. [5] which has 
been rescaled by a factor 3 for comparison (see the main text 
for details). 


to have unit length such that the lenght of the Brillouin 
zone is 2i r. For the tight-binding disks we used, however, 
the standard length a/3 between neighboring unit cells. 
As RSQA depends quadratic on momentum it needs to 
be rescaled by a factor 3. 


The only difference is that the yellow highlighted path 
A in Fig.[l3] is replaced by a subset A n of A. Ai only 
contains the most outer site of A, A 2 the two most outer 
sites and so on. The results are shown in Fig.[l7| We 
find that for the rather small disk sizes up to R = 15 
(z.e., 1350 lattice sites) we need to go at least to n = 4 
in order to get qualitative agreement with the exact 
result. Very good agreement is obtained for n = 6, 
K,( 6 \E\,E 2 ) = /C(E 1 ,E 2 ). Note that a disk with size 
R = 15 is still a very small sample. Typical samples sizes 
in experiments are magnitudes larger. This means that 
in real experiments it might be sufficient to measure the 
outermost two or three sites to obtain perfect agreement 
with the full RSQA. This opens the perspective towards 
a direct spectroscopic measurement of RSQA in a realis¬ 
tic experiments. Details of this spectroscopy aspect will 
be worked out in the future. 


Of course, the considered tight-binding disks with 
hexagonal shape are a “discrete analogue” of a disk with 
circular shape and it is therefore not surprising to find 
such good agreement. Therefore we consider also disks 
with non-circular shape as shown in Fig. [12] With the 
parameters used in Fig. 12 we find the RSQA /C(Ei,E 2 ) 
which looks qualitatively identical to Fig.[l5] With this 
final test, we have established that the RSQA will be 
present in any generic helical liquid independent of the 
considered geometry. 


VIII. CONCLUSION 

We have analyzed the spin texture of generic helical 
liquids, the gapless helical edge states of two-dimensional 
topological insulators with broken axial spin symmetry. 
Generic helical liquids feature a momentum-dependent 
rotation of the spin-quantization axis which eventually 
is responsible for the leading finite-temperature correc¬ 
tion to the otherwise quantized Hall conductance!- 18 Here 
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we considered two different topological insulator models, 
the Bernevig-Hughes-Zhang and the Kane-Mele model, 
and computed the rotation of the spin-quantization axis 
for these models in detail. For systems with translation 
invariance, we considered parameter and sample size de¬ 
pendencies. Moreover, we considered disks which do not 
possess translational symmetry anymore but rotational 
symmetry instead. We showed that in such disks the 
spin-quantization axis rotates as a function of total an¬ 
gular momentum. Finally, we showed that the rotation of 
the spin-quantization axis remains accessible even when 
neither continuous translational nor rotational symme¬ 
tries are preserved. This proves the ubiquity of the rota¬ 
tion of spin-quantization axis, which is independent of the 
considered geometry. The study of tight-binding disks 
or flakes also revealed that the information needed can 


be extracted almost locally, opening the path towards 
spectroscopic detection of the spin texture of helical edge 
states. 
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FIG. 17. Recovering the RSQA JC(Ei, E?) site-wise by “measuring” the wave functions only on the most outer site, the two most 
outer sites, etc. of a tight-binding disk, (a) E 2 ), (b) JC^(Ei, £’ 2 ), (c) JC^ 3 \Ei, E 2 ), (d) E 2 ), (e) E 2 ), 

and (f) JC^ 6 \Ei, E 2 ) are shown. For the four most outer sites (d) the well-known RSQA structure becomes visible, while it 
is fully recovered for panel (f). Parameters used: t — 1, Aso = 0.2 and Ar,i = 0.05, Ar ,2 = 0, and with disk size R — 15 
corresponding to 1350 lattice sites. 






















